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1 INTRODUCTION 


Phase Transformations in Binary Colloidal Monolayers^ 

Ye Yang,*'* Lin Fu,^* Catherine Marcoux,*^* Joshua E. S. Socolar,*^ Patrick Charbonneau,^’*^ and 
Benjamin B. Yellen"’'^’'' 


Phase transformations can be difficult to characterize at the microscopic level due to the inability to directly observe individual 
atomic motions. Model colloidal systems, by contrast, permit the direct observation of individual particle dynamics and of collec¬ 
tive rearrangements, which allows for real-space characterization of phase transitions. Here, we study a quasi-two-dimensional, 
binary colloidal alloy that exhibits liquid-solid and solid-solid phase transitions, focusing on the kinetics of a diffusionless 
transformation between two crystal phases. Experiments are conducted on a monolayer of magnetic and nonmagnetic spheres 
suspended in a thin layer of ferrofluid and exposed to a tunable magnetic field. A theoretical model of hard spheres with point 
dipoles at their centers is used to guide the choice of experimental parameters and characterize the underlying materials physics. 
When the applied field is normal to the fluid layer, a checkerboard crystal forms; when the angle between the field and the normal 
is sufficiently large, a striped crystal assembles. As the field is slowly tilted away from the normal, we find that the transforma¬ 
tion pathway between the two phases depends strongly on crystal orientation, field strength, and degree of confinement of the 
monolayer. In some cases, the pathway occurs by smooth magnetostrictive shear, while in others it involves the sudden formation 
of martensitic plates. 


1 Introduction 

Diffusionless transformations are a class of solid-solid phase 
transitions in which the crystal unit cell changes shape and 
internal structure, while keeping its stoichiometry constant. 
These include magnetostriction, where a crystal undergoes 
continuous shear, as well as martensitic transformations, 
which involve more complex particle rearrangements. Be¬ 
cause these transformations do not require long-range diffu¬ 
sion, they are fast and repeatable, and thus have been exploited 
in a number of engineering applications, including actuation in 
shape-memory alloys, hardening in steel, ^ and heat transfer 
in giant magnetocaloric materials. Some martensitic tran¬ 
sitions are irreversible, as in the quench hardening of steel, ^ 
while others are reversible, as in shape memory alloys and 
giant magnetocaloric materials. The ability to optimize dif¬ 
ferent aspects of these transformations, such as the response 
intensity and the repeatability over many cycles, requires de¬ 
tailed knowledge of how the crystal structure changes due to 
temperature and mechanical, electrical, or magnetic stresses.^ 
Because of their important technological applications, marten¬ 
sitic transformations have been extensively studied. Their 
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mesoscopic features have been probed by calorimetry, ^ acous¬ 
tic emission,^ and in situ microscopy;their microscopic 
dynamics have also been studied through numerical simula¬ 
tion of simple models. However, no experimental technique 
has yet accessed the transformation dynamics with atomic res¬ 
olution. 

Though significantly larger than atoms, colloidal particles 
represent an accessible model for studying phase transitions 
in condensed matter. Being individually resolvable by optical 
microscopy, yet sufficiently small to form thermally equili¬ 
brated phases at room temperature, micron scale colloidal par¬ 
ticles have yielded important insights into mechanical crystal 
growth, melting in confined geometries, and solid- 
solid transitions. Because of their experimental simplic¬ 
ity, monocomponent colloidal systems have garnered the most 
attention; however colloidal alloys offer a greater diversity 
of equilibrium phases and transformation pathways, though 
they are more difficult to study. 

Here, we study both fiuid-solid and solid-solid transforma¬ 
tions that occur in a mixture of magnetic and nonmagnetic 
colloidal particles immersed in a thin aqueous solution of 
magnetic nanoparticles (ferrofiuid) and exposed to an exter¬ 
nal magnetic field. This system behaves similarly to a mixture 
of point dipoles with opposite orientations, where the relative 
magnitude of the dipole moments can be tuned by changing 
the ferrofiuid concentration. We study colloidal monolayers 
by setting the thickness of the fiuid layer to be nearly equal to 
the particle diameter, thus enabling the investigation of trans¬ 
formation pathways with high spatial and temporal precision. 

A rich variety of equilibrium crystal structures can be ob¬ 
tained by adjusting the ferrofiuid concentration, the relative 
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Fig. 1 Phases of binary colloidal monolayers, (a) Magnetic particles (red) and nonmagnetic particles (blue) are depicted in the checkerboard 
crystal phase at low tilt angles and striped crystal phase at high tilt angles. Experimental images corresponding to the cartoons in (a) are 
provided in (b) the checkerboard crystal in a tilt angle of 0 = 0° and (c) the striped crystal in a tilt angle of 0 = 60°. Magnetic particles have 
average diameter of 2.8 jUm, and nonmagnetic particles have diameter of 3.1 fim. Particles are immersed a ferrofluid with 1% volume fraction 
of magnetic nanoparticles. Scale bars are 10^um. (d) Calculated minimal potential energy configuration as a function of tilt angle and negative 
susceptibility ratio, —Xm^m/Xn^n^ for a perfectly confined 2D system with equal bead size and particle number densities. Note the natural log 
scale on the horizontal axis. Gravity and image dipoles were not included in this calculation. The susceptibility ratio can be tuned by changing 
the nanoparticle concentration in the ferrofiuid. Experiments and simulations in this work were performed at —Xm^m/Xrv^rv = 1 (solid gray 
line). Note that for 3 < Xm^m/Xn^Ti ^3.7 (e.g., dashed gray line), decreasing the tilt angle transforms an equimolar striped crystal into a 
hexagonal crystal with a 2:1 ratio of magnetic to nonmagnetic particles through an intermediate chain phase. This change in crystal 
stoichiometry suggests that the solid-solid transition should be accompanied by long-range particle diffusion. 


particle density (^m/^n), and the direction of an external field. 
The minimal energy (zero-temperature) phase diagram com¬ 
puted for 65 different structures (see Fig. Id and Sec. 3.1 
for calculation details) illustrates the wealth of phases formed 
in an equimolar alloy system (^m/^n = !)• Various solid- 
solid transformations are achievable, including some diffu¬ 
sionless ones. The experiments reported here roughly follow 
the solid vertical line of Fig. Id, in which the particle moments 
are equal and opposite, and where the checkerboard crystal 
(Fig. lb) transforms to a striped crystal phase (Fig. Ic) as the 
direction of the magnetic field is tilted away from the z- direc¬ 
tion. 

Interestingly, our experimental studies of these transforma¬ 
tion pathways indicate the presence of intermediate states that 
cannot be explained by minimum energy configurations in 
perfectly confined two-dimensional (2D) monolayers. Our ex¬ 
periments, Monte Carlo simulations, and analytical calcula¬ 
tions, further reveal that the specific pathways are highly sen¬ 
sitive to the crystal orientation and to the degree of monolayer 


buckling allowed, i.e., the distance of separation between the 
plates that confine the sample. Some of these pathways ap¬ 
pear smooth and continuous over a broad range of tilt angles, 
whereas others proceed via sequences of discrete collective 
motions. 


The rest of this article is organized as follows. Section 2, de¬ 
scribes the important physical features and parameters of the 
binary colloid system. Sections 3 and 4 detail the numerical 
and experimental methods, respectively. Section 5 presents 
results from experiments and simulations, including the phase 
diagram and crystal growth, magnetos trie tive response to ap¬ 
plied fields for small tilt angles, and the martensitic phase 
transformation that occurs at large tilt angles. We briefly con¬ 
clude with remarks about equilibration times in our experi¬ 
mental system and future research directions. 


2 

















2 THEORETICAL MODEL 


2.2 Self-consistent magnetic moments and image dipoles 


2 Theoretical model 


leading to the effective dipole moment 


The magnetic particles used in experiments consist of mag¬ 
netic nanograins uniformly distributed inside an inert, spher¬ 
ical micron-sized polymer matrix. They thus behave as a 
homogeneous continuum on the micron scale, allowing for 
the assignment of an effective magnetic permeability, or 
/in, for magnetic or nonmagnetic particles, respectively. The 
nonmagnetic particles have similar composition, but do not 
contain magnetic nanograins. Both types of particles are 
immersed in an aqueous ferrofluid of magnetic nanograins, 
whose concentration tunes the average magnetic permeability 
/if. In the following, the ferrofluid is assumed to be homoge¬ 
neous. 

2.1 Dipolar potential energy 

From classical magnetostatics we know that a homogeneous 
sphere of magnetically susceptible material produces the held 
of a point dipole when exposed to a uniform magnetic held. 
The potential energy of the binary colloid system placed in an 
external magnetic held is thus modeled by taking the parti¬ 
cles to be induced point dipoles at the centers of hard spheres, 
with the effective susceptibility determined by the permeabil¬ 
ity difference between a sphere and the surrounding fluid. The 
system potential energy U is thus taken to be the sum of hard 
sphere interactions and magnetic interactions between the in¬ 
duced point dipoles. The gravitational energy associated with 
out-of-plane buckling is also considered. For N particles, we 
thus have 


N 

U = L[t/Hs(ro) + Uiiinj) + Ugin)] (i) 

i<j 

where Vij is the distance between particles i and 7 , and hard 
sphere exclusion is complete up to the particle surface at {ot + 
Oj)/! for particles of diameter cr. Dipolar interactions are 
given by the classical expression, 

Mf 

t^dd(r, 7 ) = --—^ [3 (m,- • f,y)(m^- • f^) -m,- • m^] , (2) 

4;rr^-y 

where m/ and mj are the effective dipole moments and r is 
a unit vector. The ferrofluid magnetic permeability Pf is as¬ 
sumed to be a function of the material bulk magnetic suscep¬ 
tibility Xb and the volume fraction of nanoparticles (p, 

Atf = a<o(1 + <PZb), (3) 


where po is the vacuum permeability. 

The effective magnetic susceptibility of a particle sub¬ 
mersed in the ferrofluid is given by 




f A 

VM, + 2 Mf/ 


(4) 


m; = XiViUi, (5) 

where vt is the particle volume. When pf lies between p^ ^ po 
and Pm > Po, the nonmagnetic particles are effectively dia¬ 
magnetic (Xn < 0 ), while the magnetic particles are paramag- 
netic (Xm > 0 ). 

No sufficiently accurate, direct measurements of the differ¬ 
ence in susceptibilities between each particle type and the fer¬ 
rofluid are available. We can, however, estimate the ratio of 
the two moments in an equimolar mixture by noting that the 
observed checkerboard phase is a potential energy minimum 
when the susceptibility ratio Xm^mlXn^n = — 1 and the tilt an¬ 
gle 6 is small (see Fig. Id and Sec. 3.1). In experiments, we 
focus on systems where the largest checkerboard crystals as¬ 
semble, and we thus reasonably assume that Xm^mlXn^n = — 1 
for the rest of the experimental and numerical analyses. 


2.2 Self-consistent magnetic moments and image dipoles 

The magnetic moment of a given particle should be calculated 
from the total held at its center, which is the sum of the applied 
held and the held created by neighboring particles: 


m,- = XAi 




I 


3{mj-rij)rij-m 

Anr^j 


( 6 ) 


where {H^) is the vertical (in-plane) component of the ex¬ 
ternal held in air and denotes the set of neighbors of 
particle i within a cutoff radius ^. 

Such computations are possible, but numerically expensive, 
so we first evaluate their relevance. Self-consistent calcula¬ 
tions of the effect of nearest-neighbor fields on the magnetic 
moments reveal that they generate only small contributions to 
the potential energy. The only structure with a non-negligible 
change in potential energy is the incommensurate stripe phase, 
but even in this case the zero-temperature phase boundary 
shifts by only a few degrees, from 41° to 49° (see Fig. 2). For 
computational efficiency, we therefore ignore the held created 
by other particles when calculating the magnetic moments in 
Monte Carlo simulations. 

The presence of a magnetic permeability mismatch at the 
fluid-glass interface gives rise to image dipoles. The dif¬ 
ference in magnetic permeabilities of the ferrofluid and the 
confining glass results in an additional held felt by the par¬ 
ticles. Consider a point dipole with magnetic moment m = 
(mjc, m^) located at r = (vo,yo 5 ^o) within the ferrofluid. Let 
the bottom glass slide be in the plane z = 0 and the coverslip be 
in the plane z = h. The held within the ferrofluid at r = (v,y, z) 
is thus a sum of the magnetic held of the real dipole and the 
fields of two image dipoles located at = {xo,yo,2h — zo) 
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3.1 T = 0 phase diagram 


3 NUMERICAL METHODS 


and = (xq, jo, ~^o) with magnetic moment: 

^ ] {mx, my, -mj . (7) 

\M/ + Mo/ 

The inclusion of image dipoles in potential energy cal¬ 
culations produces slightly different minimum energy phase 
boundaries, but again the phase transition sequence is not 
qualitatively affected (Fig. 2), which justifies ignoring their 
effect in simulations. 


2.3 Gravitational effects 

A gravitational contribution must be included when the exper¬ 
imental cell height is larger than one particle diameter. For a 
density mismatch pi between particles of type i and the fer- 
rofiuid, we take 

l/g(r,) = PiVigrrz, (8) 

where vi is the volume of a bead of type i and the z = 0 plane 
is at a distance CTm/2 above the bottom glass slide. In our 
experiments, the density of the nonmagnetic beads is closely 
matched with that of the fiuid, but the magnetic beads are more 
dense, hence: Pn = 0 and Pm = 350kg/m^. 


2.4 Reduced temperature 

In order to directly compare simulations, in which the temper¬ 
ature is varied, with experiments, in which the applied field 
is controlled, we define a reduced temperature as the ratio of 
thermal energy to dipolar potential energy in the system: 




(9) 


structures of equal-sized magnetic and nonmagnetic parti¬ 
cles, i.e., <ym = On- This set includes rings, chains and crys¬ 
tals that were either presented in Ref. 30, manually con¬ 
structed, or predicted by a genetic algorithm. Potential en¬ 
ergy and magnetic moment calculations were performed us¬ 
ing the same methods as in Ref. 30. Magnetic moments were 
self-consistently determined by solving Eq. (6) using a cutoff 
radius of ^ = 1.1 cr. The overall potential energy calculation 
used no cutoff radius for the rings and cutoff radii of 500(7 
and 100(7 for the chains and crystals, respectively. Gravity 
and image dipoles were neglected, and, for non-stoichiometric 
crystals, additional particles were assigned zero potential en¬ 
ergy, as in the ideal gas limit. Other parameters used in the 
calculations are presented in Table 1. Note that performing 
the energy calculations with different values of (7, |H|, and Xb 
does not change the location of the phase boundaries in the tilt 
angle-susceptibility ratio phase diagram. 
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Table 1 Parameters used in potential energy ealculations for Fig. Id 
and Fig. 2. 


where T = 298 K is the ambient temperature, is the Boltz¬ 
mann constant, (7 = ((Tm + (7n)/2 is the average particle diame¬ 
ter, and a is an effective susceptibility that accounts for exper¬ 
imental features that are not modeled directly and corrections 
to the dipole moment approximations mentioned above. 

3 Numerical Methods 

Different methods were used to extract quantitative and qual¬ 
itative information from the variants of the theoretical model. 
The numerical and simulation details are provided in this sec¬ 
tion. 

3.1 r = 0 phase diagram 

The minimal potential energy structures as a function of 6 and 
—Xm^m/Xn^n Were determined from a set 65 two-dimensional 


Comparisons of potential energies with and without self- 
consistent fields and image dipoles were made for seven types 
of structures. Here again, energy calculations were performed 
in the high-field limit, in which the gravitational energy is neg¬ 
ligible. Three buckling heights were then considered, corre¬ 
sponding to slide to coverslip separations h = 1.17 (7n, and 

1.37 (7n. 

For calculations in which magnetic moments were self- 
consistently determined, Eq. (6) was solved using the Jacobi 
method with a cutoff radius of^= 20.1 (7n. Convergence of 
the energy to within 0.1% was obtained within three iterations. 
For all cases, the total potential energy per particle was calcu¬ 
lated using a cutoff radius of Tu = 30.1 (7n. The two cutoff radii 
were chosen so as to minimize computation time while yield¬ 
ing an error of < 0.1% in the potential energy. The values of 
the parameters used in calculations are presented in Table 1. 

The values of (p and Xb were chosen such that the magnetic 
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3 NUMERICAL METHODS 


3.1 r = 0 phase diagram 
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Fig. 2 Minimal energy structures at various field tilt angles 9. Energies were calculated using different values of glass slide to coverslip 
distance h, and h= 1.00 corresponds to a perfect 2D confinement. For each buckling height, the effects of image dipoles and self-consistently 
determined magnetic moments were included in the energy calculation. In panels (a), (d), and (g) image dipoles were included and 
self-consistently determined moments were used. In panels (b), (e), and (h) no image dipoles were included and self-consistently determined 
moments were used. In panels (c), (f), and (i) no image dipoles were included and the magnetic moments were fixed by the external field 
alone. 
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3.4 Phase diagram 


3 NUMERICAL METHODS 


moments of the magnetic and nonmagnetic particles are equal 
and opposite in the limit of infinite dilution. This results in the 
checkerboard crystal being stable at 0 = 0 (Fig. Id). Note that 
because both the magnetic moments and the dipole energy are 
functions of the product (pXB , calculations performed with any 
values of (p and Xb that satisfy the magnetic moment condition 
above result in identical system energies. 

The value of jlm in experiment is not known precisely. For 
the structures observed experimentally (depicted in Fig. 2), en¬ 
ergy calculations were performed with /im = 1.5, which ac¬ 
counts well for the tilt angles at which transitions should oc¬ 
cur. 


3.2 Finite T simulations 

The phase diagram and hysteresis loop calculations were ob¬ 
tained from systems in perfect 2D confinement. For simplic¬ 
ity and efficiency, two additional approximations were made. 
First, the magnetic and nonmagnetic particles were chosen 
to have the same diameter, cr. Second, the magnetic mo¬ 
ments were not calculated self-consistently, but kept fixed (see 
Sec. 2.2 for a detailed discussion). The first approximation 
was only made for the phase diagram calculation and the sec¬ 
ond one was made for all simulations. It was empirically ob¬ 
served from comparing the results with experiments that the 
errors introduced by these approximations can be accounted 
for by rescaling the fitting factor a. 

In order to keep the error on the numerical energy calcu¬ 
lations within 0.1%, the cutoff radius for the pair interaction 
Tsim was chosen to be half the simulation box size, i.e., a given 
particle is considered to interact with ^ Nn/A particles. To 
further improve the computational efficiency, energy calcula¬ 
tions for the solid phase used Ewald summation (see, e.g. Ref. 
32, Ch. 6, App. F). For the fiuid phase, however, in the sys¬ 
tem size regime considered the approach is less efficient than 
real-space radial truncation. 

3.3 Equations of state 

The fiuid and the crystal equations of state (EOS) were deter¬ 
mined by measuring density p as a function of pressure P from 
constant NPT Monte Carlo (MC) simulations. Simulations 
with N = 400 particles were averaged over 4x10^ MC cycles 
following an initial equilibration of 10^ MC cycles. Each MC 
cycle consists of N particle moves and one volume move, on 
average. In order to efficiently treat strongly associated parti¬ 
cles at low densities (p < 10“^) and temperatures (T < 0.14), 
aggregation-volume-bias MC (AVBMC)^^ and cluster volume 
moves are substituted for half of the standard local displace¬ 
ments and volume moves. We used tavbmc = l-^CTiy as the 
cutoff for the AVBMC and r^u = TOScr^y as the cutoff for the 
cluster volume move, where Gij = (Oi + cry)/2. The maximal 


step size Ax and AlnV for displacements and volume moves, 
respectively, were tuned every 1000 MC cycles, in order to 
keep the acceptance rates of both types of moves between 30% 
and 40%. 

3.4 Phase diagram 

The system free energy was obtained through thermodynamic 
integration (see, e.g. Ref. 35, Ch. 10). Eor the fiuid, an 
equimolar mixture of ideal gases is used as the reference free 
energy 


F''^{p)=NkBT 



ln(2;rAr)1 

L V 2 / 

2N 


( 10 ) 


where the thermal de Broglie wavelength A is here set equal 
to (7. Eor the solid, an Einstein crystal with area fraction rj = 
0.739 was used instead. We found that the Einstein crystal 
limit is recovered for a spring constant of 2000^^7 / cr^ for the 
different state points studied. 

The fiuid-solid coexistence densities at a given T were de¬ 
termined from a parametric curve of the chemical potential, 
F/A + P/p, and P for each phase as a function of den¬ 
sity. Phase coexistence takes place at the intersection of the 
two curves. Note that the fiuid-checkerboard crystal coexis¬ 
tence region obtained here is consistent with an earlier, lower- 
precision calculation. 


3.5 Dynamical Monte Carlo 

The simulation results for both the hysteresis loop and the 
martensitic transformation, dynamical Monte Carlo simula¬ 
tions, i.e., using only local single-particle displacements, were 
performed at constant NVT and for N = 400. 

The hysteresis loop of Eig. 4 was determined at a fixed area 
fraction 7] = 0.46 in a 2D system. The order parameter was av¬ 
eraged over 2000 MC cycles for each T (H). The reduced tem¬ 
perature was changed by AT = 0.004 every 2000 MC cycles 
between T* = 0.023 — 0.075, AT = 0.015 every 2000 MC cy¬ 
cles between T* = 0.075 — 0.135, and AT = 0.03 every 2000 
MC cycles between T* = 0.135 — 0.375. 

Eor the martensitic transformation simulation, the full 
Hamiltonian for particles with diameter ratio CTm • CTn = 2.8 : 
3.1 (Section 2) was used. Both perfectly 2D and quasi-2D 
systems were considered. In the former, particle centers were 
constrained to only move in the z = (7n/2 plane; in the latter, 
the particles were allowed to move between hard walls sepa¬ 
rated by a distance h= l.llcTn. The external field was tilted 
from 0° to 50° with a tilt rate of 2.5°/2000 MC cycles. To 
straightforwardly compare the simulation results with the ex¬ 
periments, images and movies were colored using the same 
protocol as used in experiments (See Section 4.4 and Movie 
S4). 
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4 EXPERIMENTAL METHODS 


4.3 Martensitic transformation 


4 Experimental Methods 

Colloidal monolayers were formed by placing an equimo¬ 
lar (1:1) binary mixture of spherical 2.8/im diameter mag¬ 
netic (M-270 Dynabeads®, Life Technologies™) and 3.1 jim 
nonmagnetic particles (Fluro-Max R0300, Thermo Fischer) 
immersed in a ferrofluid (EMG705, Ferrotec, Bedford, NH) 
between a glass slide and a coverslip. Both magnetic and 
nonmagnetic particles have a size dispersity of less than 
3%. To reduce adhesion between particles and surfaces, the 
glass slides and coverslips were coated with polyethylenegly¬ 
col (PEG) by silanization with lOkD silane-polyoxyethylene- 
carboxylic acid (PG2-CASL-10k, NANOCS, New York, NY). 
The ferrofluid susceptibility was controlled by adjusting the 
volume fraction of magnetic nanoparticles in the flnal suspen¬ 
sion to ^ 1%. 

A 1.9/iL aliquot of fluid mixture was sealed in between 
the glass surfaces with Loctite marine epoxy. Uniform mag¬ 
netic flelds were applied by passing currents through air- 
core solenoids (Fisher Scientiflc, Pittsburgh, PA). Microscopy 
was performed with a DM LM fluorescent microscope (LE- 
ICA, Bannockburn, IL) using a 40X air-immersion objective 
and a combination of brightfleld and fluorescent Alter cubes 
(Chroma Technology). Videos were recorded with Retiga 
2000R camera (Qimaging, Surrey, Canada). The held ramp 
was produced by a combination of vertical and horizontal 
solenoids controlled with Lab view (National Instruments, Ver¬ 
sion 2010, Austin, TX). Image processing was performed 
using MATLAB (Mathworks, Version 2011, Natick, MA). 
Details about the image processing and data analysis meth¬ 
ods ^3,34 provided in Section 4.4. 

4.1 Phase diagram 

The fluid-solid transition was studied in a vertical magnetic 
fleld, with samples prepared at various rj. (See Fig. 5.) Each 
cooling-heating cycle lasted 512 minutes. The magnetic fleld 
was adjusted with the square root of time, which results in the 
inverse reduced temperature changing linearly with time. The 
time-dependence of the magnetic fleld over a single cycle was 
thus 

0.5v^z, 0 < r < 256 

0.5V512-rz, 256>r<512, 

where the fleld strength is in Oe and time is in minutes. An 
image of the system was taken every minute. An example of 
the fluid-solid transition in an effective cooling-heating cycle 
at T] = 0.43 is demonstrated in Movie SI. 

4.2 Magnetostriction 

Magnetostriction experiments proceed by adding a weak, in¬ 
plane magnetic fleld to a vertical fleld of 9.5 Oe, hence the 


total external magnetic fleld was 

Ho = HzZ-\-HxSm{27rft)x, (12) 

where = 9.5 Oe, = 2.1 Oe. The cycling frequency / = 
5.0 X 10“^ Hz was chosen such that the crystal had sufficient 
time to equilibrate, and it corresponds to the polar field tilt 
angle reaching 12.5° from 0° in 2.00 x 10^ s. 

The magnetostriction dependence on the polar tilt angle 
(Eig. 6d) was determined from a square section of the system 
consisting of 6 x 6 nonmagnetic particles within a crystallite 
(Eig. 6b). In order to reduce the influence of other factors, 
such as crystal size and orientation, we only analyzed crystals 
containing more than 250 particles and with their 10 axis par¬ 
allel to the v-axis (azimuthal angle Xj/ = 0°). During each of 
the eight replicates, an image of the system was taken every 
50 s, and their analysis gave consistent results. 

The magnetostriction dependence on the crystal orienta¬ 
tion, or azimuthal angle y/, was determined from crystallites 
with more than 250 particles. Crystals with similar orienta¬ 
tion (Ai/r < 5°) were grouped together and crystal growth was 
repeated until every group had at least five data sets. Within 
each group, the average relative extension ratio at 6 = 12.5° 
was calculated (Eig. 6e) from images of the system taken ev¬ 
ery 20 s. An example of magnetostriction dependence on both 
6 and y/ is presented in Movie S3. 

4.3 Martensitic transformation 

More than 10 martensitic transformation experiments were 
performed to verify the consistency of the different marten¬ 
sitic transformation pathways, and an image was taken every 
30 s for analysis in each experiment. The image processing 
and particle identification protocol described in Section 4.4 
was used to determine the particle coordinates. The particle 
bonds were then characterized using the method described in 
Section 4.5. In Eigure 7 and Movies S8, SIO, S12, and S13, 
the particles are artificially colored according to their bond 
type for visualizing the details in the martensitic transforma¬ 
tion process. 

4.4 Modified square order parameter 

Experimental measurements of d >4 were obtained from images 
taken during the cooling-heating cycle that were preprocessed 
using background subtraction and noise reduction, following 
the method developed by Crocker et al. A two-step particle 
identification was then performed: (i) nonmagnetic particles 
were identified by finding the local intensity maximum of the 
particle center; (ii) magnetic particles were identified by first 
locating the dark rim of the particle edge, and then checking 
that the intensity of the tentative particle center is within a 
properly chosen range. Visual inspection of multiple images 
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4.6 Reciprocal lattices 


5 RESULTS 


(a) (b) (c) 



Fig. 3 Illustration of the order parameter calculation, (a), Particles 
within 1.3 (7m from the particle of interest are considered to be 
neighbors and are colored gray; particles beyond that distance are 
colored white. An example experimental configuration (b), gives a 
corresponding order parameter analysis (c). Particles are colored 
according to their individual values of O 4 (Section 4.4), using the 
color scheme shown in Figure 4b. 

reveals that this identification protocol correctly determines 
> 98% of the particle coordinates. For Movies S 8 , SIO, S12, 
and S13, the remaining particles were manually corrected. 

Neighbors, i.e., particles within 1.3(7m of each other, are 
then identified. This definition deliberately excludes second- 
nearest neighbors in a square lattice. Given a particle k with 
n neighbors, the local order around particle k is analyzed in 
terms of a modified square order parameter 


^ 

'0, 

1 

n 

Y, exp(4/ej_^) 


< 

jedkW 


where denotes the set neighbors (within = 1.3 dm of 
particle k), and Oj^k is the acute angle between the x axis and 
the bond between particle k and its 7 * neighbor. This defini¬ 
tion allows us to differentiate between the checkerboard crys¬ 
tal phase (^ = 4 in the bulk, ^ = 3 at the boundary, hence 
^ 4 ,k 7 ^ 0 ) and the chain phase (^ = 2 in the bulk, ^ = 1 at the 
boundary, hence = 0 )- An illustration of this calculation 
is given in Fig. 3. The global order parameter is then defined 
as the average of < 1>4 over all particles, 

1 ^ 

= (14) 

7=1 

where N is the total number of particles in the field of view. 

4.5 Particle bond characterization 

Types of particle bonds, rather than a local order parameter, 
were used to quantify the martensitic transformation. When 
Vij < 1. 2 ( 7 / 7 , particles are deemed bonded. If the bond 

formed by a pair of like particles lies within 20 ° from the in¬ 
plane field component, then it is considered to be a martensitic 


bond; if a bond is between a pair of unlike particles, then it is 
considered to be an austenitic bond. The number and distribu¬ 
tion of bonds is used to characterize the formation of phases 
in the martensitic transformation process. 

4.6 Reciprocal lattices 

The reciprocal lattices were obtained by computing the 
Fourier transform of experimental and simulation images. To 
eliminate artifacts created by the image boundary, the real- 
space images were first processed with a Gaussian kernel 
function exp{ —[(x —+ (y — yo)^]/ 2 cr^}, where (xo,yo) is 
the center of the image and a = 1000. A Fast Fourier Trans¬ 
form (FFT) was then performed on the processed images using 
Matlab (R2011b). 

5 Results 

Results for assembly and phase transformation of the experi¬ 
mental and the theoretical models are presented in this section. 

5.1 Phase diagram in a vertical field 

To study the transformation between bulk solid phases with¬ 
out interference from edge effects and defect-mediated transi¬ 
tions, we first developed a protocol for building large, well 
ordered single crystals. This endeavor requires a quantita¬ 
tive understanding of the fluid-solid phase transition, which 
is achieved by directly comparing experimental results with 
the phase diagrams obtained by Monte Carlo simulations. 

Figure 4b shows the simulated phase diagram in the plane 
of reduced temperature T* and area fraction rj for dipolar 
hard spheres under 2D confinement. The corresponding ex¬ 
perimental results are obtained by continuously varying the 
magnitude H of 3. vertical magnetic field for samples with dif¬ 
ferent 7 ]. A bond-orientational order parameter ^4 (defined in 
Section 4.4) tracks the first-order freezing and melting transi¬ 
tions as the field strength is cycled between 0 and 8 Oe (Fig. 4a 
and Movies SI and S2). The center of the hysteresis loop is 
used to estimate the melting point at a given T]. In agreement 
with simulations, the melting field strength is roughly inde¬ 
pendent of T] at 0.1 < 7] < 0.4, as evidenced by the plateau of 
Fig. 4b, but decreases at high rj. Matching the experimentally 
observed transitions to the computed phase diagram at several 
different area fractions gives a fitting factor a = 2.4 ± 0.3, for 
Eq. (9) 

In the fluid-solid coexistence regime, we expect that large 
regular crystals should form at values of H and rj near the neck 
of the phase diagram around T* =0.12. The small density 
difference between the two phases in that regime helps anneal 
crystal defects. Preparing monolayers with 7] >0.60 was too 
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5.2 Tilting the held at small angles: magnetostriction 


(a) 



Area Fraction of Particles,^ 


(a) (b) (c) 





77=20% 

r]=33% 

77 = 41 % 

i) 


(f) 



m 

H=2.5 0e 

H=5 0e 

H=7.5 Oe 


Fig. 5 Crystal growth in a constant vertical field. The dependence 
of crystal growth on area fraction is demonstrated for (a) rj = 0.20, 
(b) rj = 0.33, and (c) T] = 0.41 at a field strength of H = 6 Oe. The 
dependence of crystal growth on field strength is demonstrated for 
(d) H = 2.5 Oe, (e) H = 5 Oe, and (f) // = 7.5 Oe at a particle area 
fraction of T] ^ 0.51. In each experiment, the magnetic field was 
held constant for six hours. The largest crystals were found to grow 
at T] ^ 0.51 and H = 5 Oe. 


Fig. 4 Fluid-solid phase transitions of binary colloidal monolayers 
in a vertical field, (a) Experimental and simulated hysteresis curves. 
Up (down) pointing triangles correspond to experimental data, while 
the red (black) dashed curve corresponds to the numerical 
simulation results for the increasing (decreasing) magnetic field. 

The triangle colors represent the <I >4 bond orientation order 
parameters, following the color scheme used in (b). Simulation 
results (dots connected with dashed lines) are superimposed, (insets) 
Typical experimental images with each particle colored according to 
its individual <I >4 value. Hysteresis curves are obtained with a field 
ramp rate that grows with the square-root of time, corresponding to 
a linear temperature ramp, (b) Phase diagram from simulations 
shows the fiuid (FI), fluid-checkerboard crystal coexistence (Fl/Cb), 
checkerboard crystal (Cb), and hexagonal striped crystal (Hex) 
phases. Black dots and connecting lines indicate the simulated phase 
boundaries. The colored bars represent the system <I >4 bond order 
parameter taken from the heating/cooling cycles of part (a) at four 
different area fractions. The left and right colored bars at each 
particle concentration represent the cooling and heating processes, 
respectively. Particles are immersed in a ferrofiuid with 1 % volume 
fraction of magnetic nanoparticles. The numerical fitting parameter 
is a = 2.4 zb 0.3. The color scheme is the same as Fig. 3. 


difficult to achieve because of the high viscosity of these so¬ 
lutions, but reasonably large crystals were nonetheless grown 
near the phase boundary at r] = 0.51 and if = 5 Oe (Fig. 5e). 
Crystals with more than 500 particles, such as those in Fig. 1, 
formed within hours and were taken as the starting points for 
studies of magnetostriction and martensitic transformations. 

5.2 Tilting the field at small angles: magnetostriction 

Tilting the magnetic field away from the vertical induces mag- 
netostrictive compression for tilt angles 0 > 0°. Experimen¬ 
tally, a strong vertical field (9.5 Oe) is applied, in order to re¬ 
duce thermal motion, and a weak in-plane sinusoidal magnetic 
field is supplemented. The in-plane component was set to os¬ 
cillate at 5.00 X 10“"^ Hz, so that the crystal had sufficient time 
to equilibrate. Repeating experiments at rates corresponding 
to oscillation periods between 2.00 x 10^s to 5.00 x 10^s gave 
comparable results (Fig. 6f), suggesting that our experiments 
were performed at (or near) thermodynamic equilibrium. 

Figure 6 illustrates the magnetostriction coefficient for dif¬ 
ferent 0 and crystal orientations, if/, relative to the in-plane 
component of the field. The components of the 2D strain 
tensor are given by £xx = (ai — a\)/a\, £yy = {bf — bi)/bi, 
£xy = cjai, and £yx = d/bi, where aj and b^ are the unstressed 
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Fig. 6 Magnetostriction of binary colloidal crystals. Experiments were performed in a magnetic field with a constant vertical component of 
9.5 Oe and a sinusoidally varying horizontal field with a magnitude of 2.1 Oe and a frequency of 5 x 10“^ Hz. (a) General deformation of a 
rectangle. The black rectangle with lengths ai and bi represents the unstressed state, and the blue shape represents the deformed state with 
lengths, af, bf, c, and d. An example experimental region used for measuring the magnetostriction effect is illustrated in both (b) real-space 
and (c) reciprocal-space. The Joule magnetostriction coefficient, defined as the extension ratio E is plotted in (d) as a function of the tilt angle 
for fixed crystal orientation Xj/ = 0°, and (e) as a function of the crystal orientation for a fixed tilt angle 6 = 12.5°. The insets in (d) illustrate 
an exaggerated crystal compression for small and large tilt angles for a fixed i/r = 0°. The insets in (e) illustrate the field direction relative to 
the crystal orientation for y/ = 2.5° and yr = 32.5°, respectively, with fixed 6 = 12.5°. The extension ratio as a function of tilting frequency is 
illustrated in (f), where 1// corresponds to the period of the in-plane sinusoidal wave. The black arrow indicates the tilting frequency that is 
used for obtaining the data in panels (d) and (e). 


(6 = 0°) lattice constants in the x and y directions, respec¬ 
tively, and Gf, bf, c, and d are measured in the final state 
(Fig. 6 a), using the original lattice in the experimental image. 
The non-magnetic particles at the corners of the rectangle de¬ 
picted in Fig. 6 b were used for the measurements. Both ai and 
bi are approximately equal to V2(am + crn)/2. We find that 
dilation (e^x + rotation (e^y — £yx)^ and transverse shear 
(Sxy + £yx) are all negligible. Figure 6 thus only depicts the 
extension ratio (or longitudinal shear) E = which re¬ 

sults from Joule magnetostriction. Negative values indicate 
that the crystal is compressed along the direction of the in¬ 
plane field. For y/ = 0° (a field aligned with the 10 crystal 
direction) the magnetostriction is maximal, exceeding 10 %, 
which is an order of magnitude larger than that observed in 
giant magnetostrictive materials; ^ by contrast, for y/ = 45° (a 
field along the 11 crystal direction), magnetostriction is sup¬ 
pressed (see Movie S3). To test the accuracy of our measure¬ 
ment of Joule magnetostriction, the reciprocal lattice obtained 
via Fast Fourier Transform (FFT) of the real-space images was 
also analyzed. The peaks at the comers of the rectangle in 


Fig. 6 c were used to calculate the magnetostriction coefficient. 
The difference between the real- and reciprocal-space results 
is within experimental uncertainty. 

5.3 Large tilt angles: martensitic transformations 

To explore the transformation from checkerboard to striped 
crystal, we applied a rotating magnetic field with constant an¬ 
gular frequency and field strength. For 0 > 15°, the crystal un¬ 
dergoes a diffusionless transformation through a combination 
of compression, shear, and slip along the 10 or 11 lattice di¬ 
rections. Video microscopy and numerical simulations reveal 
a variety of possible combinations of these different elements 
(Fig. 7). For a fixed T (constant H), the chosen pathway in¬ 
deed depends on the initial crystal orientation, field strength, 
and degree of planar confinement. Note that the latter two fac¬ 
tors are related because out-of-plane motion is suppressed near 
the melting transition, where the gravitational cost of buckling 
is large compared with the energy gain of forming dimers. 

To characterize the different experimental transformation 
pathways, we first consider the behavior of systems where the 
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Fig. 7 False color experimental and simulation images of the martensitic transformation from a checkerboard crystal phase to a striped crystal 
phase for different crystal orientations and field strengths. Each row shows a sequence of images from a single experiment or simulation in 
which the tilt angle of the field was monotonically increased. For visualization purposes, the particles with at least one martensitic bond are 
falsely colored red (pink) if they are magnetic (nonmagnetic), while particles with only austenitic bonds are colored dark blue (light blue) if 
they are magnetic (nonmagnetic). Simulated transformations under perfect 2D confinement and experimental transformations at both strong 
(12 Oe) and weak (5 Oe) fields are shown for each of the two crystal orientations i/r = 0° and y/ = 45°. By setting /i = 1.11(7„ in simulations 
(Movies S9 and SI 1), better agreement with experiment is obtained, which highlights the role of confinement in the transformation dynamics. 
In all experiments, the tilt rate of the field 0.008°/s. 
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motion is perfectly confined to a 2D plane, which are only 
accessible in simulation. In this case, for all field strengths 
the transition at i/r = 0 ° proceeds through continuous longi¬ 
tudinal shear (see Fig. 7(row a) and Movie S4 for the high 
field cases, and Movie S5 for the low field cases), while at 
Xj/ = 45° zigzags form and gradually straighten into horizon¬ 
tal stripes (see Fig. 7(row d) and Movie S 6 for the high field 
cases, and Movie S7 for the low field cases). These pathways 
are qualitatively different from those observed in experiments 
(see Fig. 7(rows b, c, e, and f)), in which perfect 2D confine¬ 
ment is not accessible. Although we cannot precisely quantify 
the amount of experimental buckling, we can estimate it by 
comparing the distance between adjacent particle centers rela¬ 
tive to the diameter of the particle. This analysis suggests that 
in most experiments the buckling angle between two adjacent 
particles was in the range of 5° to 15° relative to the plane. 
When a similar amount of out-of-plane motion (buckling) is 
allowed in simulation, we find markedly improved agreement 
with experiment and thus conclude that confinement plays a 
key role in tuning the transformation pathway. 

Each row in Fig. 7 shows the crystal transformations for 
different tilt angles 0 and crystal orientations y/. The first 
and fourth rows depict the simulated transformation for strict 
2D confinement. The second and fifth rows depict the exper¬ 
imental transformation in a strong field. The third and sixth 
rows depict the corresponding experimental transformations 
in a weak field. 

Experimental observations of the transition from checker¬ 
board to striped phase for two field strengths and two initial 
crystal orientations are shown in Eig. 7. The experimentally 
observed transformation pathway for a strong field and i/r ^ 0 ° 
proceeds first through magnetostrictive compression, and then 
by shear, resulting in the diagonal lines of dimers visible in 
Eig. 7b(iii). Adjacent dimer lines then align and coalesce 
to form the striped phase of Eig. 7b(iv). (See experimental 
Movie S 8 and simulation Movie S9, which allows buckling, 
unlike the simulation of Eig. 7a(iii)). These transformations 
are smooth and homogeneous. By contrast, the transforma¬ 
tion for Xj/^ 45° proceeds through a sequence of abrupt slips, 
giving rise to the vertical lines of nearly horizontal dimers vis¬ 
ible in Eigs. 7e(ii) and 7e(iii). In most cases, the dimer line 
nucleates at a defect site and is facilitated by a local lattice 
expansion, which opens up free space for adjacent columns 
to shift and for dimer lines to zip up. Isolated dimer lines 
first form throughout the crystal, typically appearing every few 
columns. Dimer lines and monomer lines then coalesce into 
trimers or quadrimers (see experimental Movie SIO and simu¬ 
lation Movie Sll), depending on the spacing between nearby 
dimer lines, and finally form horizontal stripes (Eig. 7e(iv)). 
The sudden coordinated line slips result in structures that re¬ 
semble martensitic plates. Eor weaker fields, the features 
observed at strong fields remain qualitatively discernible, but 


the coherence of the dimer lines is weakened, as displayed in 
Eigs. 7c(iii) and 7f(iii) (see experimental Movies S12 and SI3, 
and simulation Movies S14 and SI5). We attribute this effect 
to the increased importance of gravity relative to the dipole 
interactions, which suppresses buckling and thus hybridizes 
the pathways with those of perfectly confined systems. The 
perfectly 2D confined systems, which can only be probed via 
simulations, show little evidence of a dimer phase, as may be 
expected given that dimer formation in the experimental sys¬ 
tem involves buckling (see Eigs. 7a(iii) and 7d(iii)). 

Eourier transforms of the experimental images (see insets 
in Eig. 7 and Movies S 8 , SIO, S12, and S13) confirm the ro¬ 
tational symmetries of the different intermediates and show 
evidence of ordering at wavenumbers present in neither the 
initial checkerboard not the final striped phase. Eor Xj/ ^ 0°, 
the system is first continuously sheared. Let qi = (^i,0) 
and q 2 = ( 0 ,^ 2 ) be the wave vectors of the fundamental 
Eourier peaks. During the magnetostriction, qi grows and q 2 
shrinks. As dimer lines form, diffuse Eourier peaks appear 
at ±(—qi +q 2)/2 (Eig. 7b(iii) and Movie S 8 ), indicating the 
emergence of a dimer configuration with a unit cell twice as 
large as the (sheared) checkerboard. As the field is tilted fur¬ 
ther, however, both the peaks at (—qi ± q 2)/2 and at qi are 
extinguished, a signature of the striped phase. While trimer or 
quadrimer patterns may also exist, we see no clear evidence of 
ordering at the wavenumbers associated with them. 

Eor Xf/ ^ 45°, the rectangular symmetry of the Eourier pat¬ 
tern is retained throughout the transformation, and Eourier 
peaks remain fixed in position, indicating that there is no over¬ 
all shearing or compression of the lattice (see Movie SIO). 
Let q 3 = (1,1)^ and q 4 = (1,-1)^ be the wave vectors of 
the fundamental Eourier peaks in the checkerboard phase in 
this orientation. The presence of an ordered dimer config¬ 
uration is then signaled by the emergence of Eourier peaks 
at (qs + q4)/4 ± (qs - q4)/2 = (1/2, ±1)^ (see Fig. 7e(iii) 
and Movie SIO), which then vanish as the dimers join to form 
longer chains and eventually extended stripes. 

Eourier images of the striped phase display two columns of 
peaks at q^ ~ L2q. The peaks at nonzero values of q^ show no 
evidence of different spacings for the two stripe types. Though 
potential energy calculations indicate that an incommensurate 
phase is the ground state for a tilt angle of 50°, the Eourier 
images indicate that at finite temperatures the system does not 
display this feature. 

6 Conclusions 

We have observed several different types of fiuid-solid and 
solid-solid transformations in a binary system of magnetic and 
nonmagnetic spheres. By exploring the mechanism behind 
fiuid-solid transformation, we were able to optimize the exper¬ 
imental conditions so that a large single domain checkerboard 
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lattice can grow within a few hours. By probing the solid-solid 
transformations with a tilted magnetic field, we also found 
that the diffusionless transformation from a checkerboard to 
a striped phase can be selected by choosing the applied field 
strength and the orientation of the initial crystal with respect 
to the tilt direction of the applied field. For high fields, the 
different pathways both pass through intermediate states dom¬ 
inated by buckled dimers before reaching the striped phase. 
Because of the slow relaxation timescales along the transfor¬ 
mation pathway, we suspect that the experiment is not per¬ 
formed in the quasistatic regime. Dimer formation and align¬ 
ment may thus be a far from equilibrium feature. We note, 
however, that potential energy calculations suggest that at least 
two different dimer phases may be stable at certain tilt angles 
(see Fig. 2). At low temperatures (or high field strengths), 
phases of trimers, quadrimers, etc. may also play a role if the 
field can be tilted sufficiently slowly. For the relatively rapid 
tilt rates employed in our experiments, the high degree of or¬ 
der observed quite far from equilibrium is rather remarkable. 

In general, the ability to switch among different pathways 
could enable the control of functional characteristics of engi¬ 
neered materials, including their capacity for heat exchange 
and susceptibility to shape change along the transformation. 
The binary magnetic system described here also opens the way 
for studying other dynamical processes in alloys. For instance, 
changing the ferrofiuid susceptibility would enable the inves¬ 
tigation of solid-solid transitions that require long-range dif¬ 
fusion, such as the transformation from a striped crystal to a 
hexagonal crystalwith = 2 (see Fig. Id). 
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